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Abstract 

The dynamics of a one dimensional quantum walker on the lattice with two internal 
degrees of freedom, the coin states, is considered. The discrete time unitary dynamics is 
determined by the repeated action of a coin operator in U (2) on the internal degrees of 
freedom followed by a one step shift to the right or left, conditioned on the state of the 
coin. For a fixed coin operator, the dynamics is known to be ballistic. 

We prove that when the coin operator depends on the position of the walker and is 
given by a certain i.i.d. random process, the phenomenon of Anderson localization takes 
place in its dynamical form. When the coin operator depends on the time variable only and 
is determined by an i.i.d. random process, the averaged motion is known to be diffusive 
and we compute the diffusion constants for all moments of the position. 

1 Introduction 

The dynamics of Quantum Walks (QW for short) have become a popular topic in the 
Quantum Computing community as the simplest quantum generalization of classical ran- 
dom walks, see for example the reviews [3], [20], [23]. In the same way classical random 
walks play an important role in theoretical computer science, typically in search algorithms, 
QW provide a natural and fruitful extension in the study of quantum search algorithms, 
see e.g. [32], [4], [26] and the review [30] and references therein. On the other hand, QW 
also appeal to physicists interested in quantum dynamics. Indeed, QW can be considered 
as simple discrete dynamical systems governed by an effective unitary operator, not nec- 
essarily given as the exact exponential of i times a microscopic Hamiltonian. For a few 
models of this type, see e.g. [1], [27], [25], [9], [29], and [6], [8], [13], [16] for their mathemat- 
ical analysis. Moreover, there are recent experimental realizations of QW dynamics: [18] 
showed that cold atoms trapped in optical lattices exhibit a QW for suitably monitored 
optical lattices and [34] show that the same is true for ions caught in monitored Paul traps. 

While several types of QW have been defined and studied in different contexts, we will 
focus on the simplest one dimensional, discrete time QW on the lattice, defined in analogy 
with the classical random walk on the lattice. Consider a quantum walker on the lattice 
carrying two internal degrees of freedom (spin states), called the coin states in this context. 
The unit time step dynamics is defined by the action of a coin operator in U (2) on the 
internal degrees of freedom followed by a one step shift to the right or left, conditioned 
on the state of the coin degree of freedom, see (2.5) below. As is well known, when the 
coin operator is identical at each time step, the QW typically exhibits ballistic dynamics 
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due to translation invariancc. Moreover, if the coin operator is ehosen at random at each 
time step, which yields a non autonomous random dynamical system, then, typically, the 
averaged motion is diffusive, see e.g. [23], [24], see also [31]. 

In this paper we address the time-homogeneous case where the coin operator depends 
on the position of the walker on the lattice and is given by a sequence of random matrices 
in U{2). This defines a QW in a random environment. Such situations were first consid- 
ered numerically in [19], [33] for continuous time QW on graphs, that is, for an evolution 
operators generated by some random discrete Hamiltonian. The outcome of these numeri- 
cal works provides evidence that Anderson localization, in a strong dynamical form, takes 
place. Discrete time QW in random environments were then considered in [21]. By con- 
trast, the author establishes that a certain choice of random coin operators does not lead 
to Anderson localization. The random coin operators in [21] are given by a fixed unitary 
matrix whose diagonal elements carry random phases. 

We consider here a family of i.i.cl. random coin operators characterized by general 
requirements on the amplitude and transition probabilities to the right and to the left 
expressed as Assumptions (a), (b) and (c) below. This family is indexed by a real deter- 
ministic parameter and the randomness is determined by i.i.d. phases carried by all matrix 
elements of the coin operator, see (2.8). We prove that for all values of the deterministic 
parameter, dynamical localization takes place everywhere in the spectrum, for almost all 
realizations of coin operators. 

The lack of Anderson localization in the model considered in [21] is explained by the 
existence of a gauge transformation that fully eliminates the randomness of the model. The 
set of spatially random unitary matrices we consider coincides with the one considered by 
[24] in the study of temporal randomness, in one space dimension. 

For completeness, we briefly reconsider the study [24] of the non-autonomous case where 
the coin operators are random in time. W; relate the avciragcd motion to that of a persistent 
random walk and provide a simple proof of the fact that for all L € N, the moments of order 
2L, at time n, behave as D{2L)n^, for large n, with an explicit formula for the diffusion 
constants D{2L) > 0. 

2 Setup and Main Results 

The Hilbert space of pure states 

n = c^®i^iz) (2.1) 

represents two internal degrees of freedom, also called a coin, with Hilbert space C^, and 
the walker whose position Hilbert space is We fix a canonical basis of denoted 

by It); find the position basis consisting of vectors denoted by n G Z (eigenvectors 
of the position operator, the operator of multiplication by the variable n). 

The dynamics of the system is composed of discrete steps, each step consisting of a 
unitary evolution of the coin (operator C on C^) followed by the motion of the walker, 
conditioned on the state of the coin. The latter step is determined by the action 

\t)<»\n) ^ \t)^\n+l) (2.2) 
\i}<»\n) ^ |;)0|n-l) (2.3) 

extended by linearity to H. This means that if the coin is pointing up the walker will move 
to the right one step, and if the coin points down the walker moves to the left. The action 
of (2.2), (2.3) is implemented by the unitary operator 

S = ^{Pt<8\k + l)(/e| + ® |/e - l)(fc|} 
fcez 
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where we have introduced the orthogonal projections 



P-r = |t)(t| and P; = U)U|. (2.4) 

The one step dynamics consists in tossing the quantum coin and then performing the 
coin dependent shift 



U = S{C(E) I) with C = 



a b 
c d 



s.t. C* = C-\ (2.5) 



Hence, if one starts form the state 1 $5 \n), the (quantum) probabihty to reach, in one 
time step, the site |n + l) equals |ap whereas that to reach |n— 1) equals 1— |ap. Similarly, 
starting form | J,) |n), the probability to reach the site |n — 1) equals \a\^ and that to 
reach \n + 1) is 1 — |n|^. The evolution operator at time n reads U". 

Despite the similarity of this dynamics with that of a classical random walk, there is 
nothing random in the quantum dynamical system at hand. The dynamics is invariant 
under translations on the lattice Z, which hints at ballistic transport properties. 

More precisely, let X = I(E)x be the operator defined on its maximal domain in C'^(E)P{'Z), 
where x is the position operator given by x\k) = k\k), for all k G Z. For any L > and 
any ^ in the domain of X^, we define 

(X^)vt(n) := {-f, [/-"X^?7"*> . (2.6) 

The analog definition holds for (|X|^)$(n). We recall in appendix the proof that in our 
one dimensional setup we have 

Lemma 2.1 (Deterministic walk) Let ^' belong to the domain of X^. Then 

Um ^ '-^^ =B>0 

n—>oo n 

with B = iff C is off diagonal. 

A QW in a non-trivial environment is characterized by coin operators that depends on 
the position of the walker: for every k G Z we have a unitary Cfe on C^, and the one step 
dynamics is given by 

U = Y1 {^tCfe |fc + + PlCk |fc - . (2.7) 

fcez 

We consider a random environment in which the coin operator Ck is a random element of 
U{2), satisfying the following requirements: 

Assumptions: 

(a) {Ck\kei f^rc independent and identically distributed ?7(2)-valucd random variables. 

(b) The quantum amplitudes of the transitions the right and to the left are independent 
random variables. 

(c) The quantum transition probabilities between neighbouring sites are deterministic and 

independent of the site. 

The first assumption means that the sites fc G Z are independent. If the walker is local- 
ized at site n, then the state of the system is of the form <^(8) |n) for some normalized ip 
The probability amplitudes for transitions to states x (g) |n + 1) and x (g) |n — 1) (x € 
normalized) are (X) t) (t) C'n'P) and (X) 4-) (i) C'n'P)) respectively. The second requirement 
says that these transitions are statistically independent. Finally, the last requirement says 
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that randomness appears as phases (see after (2.5)) and that, in a certain sense, we remain 
close to a classical asymmetric random walk on the lattice. 



We show in Lemma 2.5 below that under assumptions (a), (b) and (c), we can consider 
without loss of generality 



"It 



with < r, t < 1 and + = 1 



(2.8) 



and with {w]J}fcez U {w^ifeez i-i-d. random variables. Thus, the action of the random 
operator Uc^ is 



(2.9) 



= e-'""t|t) |n+ 1) +6^'"-^!;) ® |n- 1) 
U^\i)(^\n) = -e-''""r|t)<»|n+l)+e-'"if|i)(^|n-l). 
The elimination of the randomness in the model [21] is a consequence of Lemma 2.5 also. 

To define the random phases properly, we introduce the probability space (f2,J^, P) 
where fl = T^, (T = R/2nZ), T is the a-algebra generated by cylinders of Borel sets, 
and P = ®k,^iii, where /i is a probability measure on T. We note expectation values with 
respect to P by E. 

The random variables J\. and o;^ on (il, P) are defined by 



cj^ : -> T, uj\,= uj2k 



A realization is thus denoted by w = (• • • ,Lu_-j^,uj_-j^,ujQ,ujQ,uj{,ujj^ 



,•••) e 



(2.10) 



Main result. Our main result is Theorem 2.2 and its corollary. Let Uoj be the one 
step dynamics of a QW in a random environment defined by (2.7) with Ck, k gZ given by 
(2.8), where {ix)*}j,gN,#e{t,4.} i.i.d. random variables defined in (2.10), distributed 

according to a measiire /i on T. 

Theorem 2.2 (Spatial disorder) Suppose that fi is an absolutely continuous measure 
diJ,{u)) = T{ui)dui, with density r e L°°{T), and that the support of n contains a non-empty 
open set. Then, for any r G (0, 1), there exist C < oo, a > such that for any j,k G Z 

and any a,T G {t, i} 



E 



sup 

/6C'(S),||/||oo<l 



{(T®3,f{Uu,)T<S)k) 



(2.11) 



Remarks: 

1) The extreme cases r = 1, t = and r = 0, t = 1 lead to deterministic results that are 
addressed in the remarks following Lemma 2.5. As seen from (2.9), when r = 0, the up 
and down components are independent and propagate respectively to the right and to the 
left. We get ballistic behaviour for any deterministic choice of phases. When t = Q, the 
walker oscillates back and forth, switching its coin state. Hence localization takes place for 
any deterministic choice of phases. 

2) Specializing the result to the function f{z) = , z G S, (2.11) implies the following 
almost sure result on the evolution in time of the quantum moments of the QW, see [16]. 

Corollary 2.3 There exists a set VIq of probability one, such that for any lo G Qq, any 
L > and for any 4* € (8) '^(Z) of finite support, there exists C^i < oo with 



sup(|X|^)*(n) < 



(2.12) 
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Remarks: 

1) Another Corollary of Theorem 2.2 is that the spectrum of Uu is pure point, almost 
surely, see [16]. 

2) The proof of this localization result leans on the fact that the unitary operator Uuj can 
be viewed as a doubly infinite five-diagonal band matrix on /^(Z), a set of operators first 
considered in [8]. Moreover, the randomness appears in such a way that it is possible to 
adapt the Aizenman-Molchanov method [2] developed for the study of Anderson localiza- 
tion in the self-adjoint case to this unitary setup, along the lines of [16]. The random 
operator at hand differs from those considered in [8, 16] which forces us to revisit the ar- 
guments based on the analysis of products of random transfer matrices and the associated 
Lyapunov exponents coupled with fractional moment estimates of finite volume restrictions 
of the associated resolvent. 

For the sake of comparison, we compute the asymptotics of the expectation of all 

(integer) moments of the position operator for a dynamics that is random in time. We 
follow [24] and consider the random evolution operator at time n given by 

U^{n,0) = UnUn-i---U2Uu where Uk = S{Ck®I), (2.13) 

where the i.i.d. random variables {Ck}ke'z defined by (2.8). The random phases are 
assumed in this case to satisfy 

E{e-''^l) = E(e-'"^) = 0, (2.14) 

but their common probability measure ji is otherwise arbitrary. We extend the results of 
[24] on the diffusive behavior of the QW in this setup by proving the following 

Proposition 2.4 (Temporal disorder) Assume the evolution operator is given by (2.13) 
with random phases (2.10) distributed according to a measure dfi such that (2.14) holds. 
Let ^0 = Vo ® \Q) be of norm one. Then, for any L e N, we have 

^^mX^hM)^ (2.15) 
n->oo n I 

where 

^^^^ = { 1 • 3 • 5 • • • (L - l)(iVr-')^/' |i 2 e^l. ^^'^^^ 

Remark: 

As observed in [24], assumption (2.14) allows one to express the averaged motion as a 
persistent classical random walk on Z, which is known to display a diffusive behaviour. 
Considering the associated generating function and analyzing its large times asymptotics, 
we are able to make explicit all diffusion constants. 

Before we turn to the proofs of these results, we briefly investigate the invariance prop- 
erties of the model stemming from its structure (2.7) and assumptions (a), (b) and (c). 



2.1 Structure of 

Lemma 2.5 Under the assumptions (a), (b) and (c), the operator defined by (2.7) is 
unitarily equivalent to the one defined by the choice 



"It 

"kr 



• t 
^tt 



where <t,r <1 and -h i = 1 



(2.17) 
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and {ijj\]kei U {io\}ke7. '""e iid random variables defined by (2.10), up to multiplication by 

a global deterministic phase. 

Proof. We first note that assumption (a) implies that we can consider one Ck at a time 
and (b) and (c) imply that the randomness appears as phases only (see after (2.5)). The 
indcpcindcncx; of the right and left probability amplitudes implies that the rows of Ck are 
independent random variables. 

As Cfe must be unitary for all realizations, the scalar product of its columns equals zero. 
This shows that the random phases of the elements on the same line are identical, i.e. 



■ 

"kc 



"lb 







a b 
c d 



(2.18) 



with C = 



a b 
c d 



unitary and independent of k. Parametrizing C as 



C 



ire 



n 



te- 



with < r, t < 1 and a, 7, 6* e T, 



(2.19) 



one sees that at the cost of multiplying by e'^, we can assume ^ = 0. Moreover, with 

S = ^ , we compute 

— le ' 



te- 



(2.20) 



Up to unitary equivalence, we can assume that C has the form of the r.h.s. of (2.20). 

To get rid of the last phase a, we make use of the representation of in the ordered 
basis {. . . , It) (8i |n - 1), li) (g) |n - 1), |t) <S) \n}, \ i) (g) |n), . . .} as the band matrix 



D^S, with S 



te- 






— r 



te 



r 


te'" 














—ia 


— r 



te- 



te'" 





(2.21) 



where (upon relabeling the indices of the random phases) D^j is diagonal with i.i.d. entries. 



diag(. 



%e 



. . .), and the diagonal of S consists of zeroes. Our convention 



is to fix a labeling of the canonical basis e„ of that space so that the odd rows contain 

r, te'" and the even rows contain te~'", — r. 

Following [8], we introduce the unitary operator V defined by Vcn = e"»"e„, where 
€ T, n e Z. We compute, for any fc G Z, 



V-'U^Ve2k 
V-^U^Ve2k 

Choosing for any fc e Z 



gi(C2i!-C2i;-l)g — il^2i;-l j,g2^_]^ 



+=ie-'"e2fe+2 (2.22) 



V~'^U^Ve2k+i = e'('^^'=+i"^^'=-i)e"'"^'=-^te'"e2fe-i - e'(^^''+i"^^'=+^)e"''^='=+=^re2fe+2. 



C2fe+2 = C2/S+I 



-ka. 



(2.23) 
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yields the result. m 
Remarks: 

1) The last argument of the proof above can be adapted to show that the randomness of 
the model of QW in random environment considered by [21] can be gauged away, which 
explains the absence of localization. Indeed, this model is characterized by random coin 
matrices of the form (with t = r = 1/ v^) 



r -te-'"" 



where {wfejfeez are i.i.d. random variables on T. (2.24) 



The corresponding operator then satisfies 



re2k+2. (2.25) 



Choosing Co = 0, C-i 



= and, for any k G Z+, 

C2fe+2 = -(W2fe+W2fe-2H hWo), 

C-2fc = i^-2k + W_2fe-2 H + W_2, 

C2fe+1 = —{i^2k + W2fc-2 H + Wo), 

C-(2fe+l) = W_2fe + ti;_2fe-2 H +'^-2> 



(2.26) 



the the deterministic coin operator C 



we see that Uu, is unitarily equivalent to the deterministic operator Uq, characterized by 

' t r ' 
r ~t 

2) We will consider from now on to be defined by (2.21) (with a — 0) acting on the 
Hilbert space P(Z), with canonical basis e2k = \ \) ® \k), e2k+i = | i) <8i |fc). We note 
that the position operator x on P(Z) is related to X defined by (2.6) on by 
{x — l)/2 < X < x/2, so that dynamical localization is equivalent in both representations. 

3) In case r = 0,t = 1, (2.21) shows that Uu: is equivalent to a direct sum of two shifts, 
which, in other words, corresponds to Cfc = I, fc € Z. This leads to ballistic behaviour. In 
case r = l,t = 0, the two-dimensional subspaces span{e2fe-i, e2fe}, fc € Z, are invariant, 
which forbids any kind of transport. 

4) The appealing five diagonal representation (2.21) of QW on N is also related so called 
CMV matrices associated with certain orthogonal polynomials on the unit circle, when 
restricted to l^(Z+). This fact is used in [10] to study certain properties of deterministic 
QW. 

Finally note that the structure of S can be described also via the two unitary operators 
("even" and "odd") 



t 



and Bo 



1 

1 



(2.27) 



S is the product of SoSg, where Sg is block-diagonal, with blocks Bg placed so that the 
(1, 1) entry of Be lies on even indices of the diagonal, and So is block-diagonal with blocks 
Bo placed so that the (1,1) entry of Bo lies on odd diagonal elements. Obviously Sg and 
So are unitary, and thus so is S and hence U^- 



3 Proofs 

We make use of the structure just described to define finite volume approximations of S 
and thus U^,. 
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3.0.1 Finite and semifinite volume truncation 

On P{{2no,2mo}) we define the matrices 



1 



Bo 



Br, 



B, 



B, 



1 



(3.1) 



where each 1 is a 1 x 1 block, and the 2x2 blocks Bo,Bf. are given in (2.27). The finite 
volume propagator is 

" r t 
r t 
t -r 

r i 

1 -r 

t -r 



r t 



1 

t -r 



= D^S, with S = SoSe 



. (3.2) 



Here, = diag(c 



Similarly we define = D^S^ acting on the Hilbert spaces ?^({2no, . . .}) and Z^({. . . , 2mo}). 
The matrix is defined as in (3.2), but where the 2x4 blocks are repeated indefinitely 
towards the bottom right. Similarly we define , and the are simply the correspond- 
ing semifinite truncations of the diagonal matrix D^.. Analogous definitions hold for the 
Hilbert spaces l'^{{2no + 1, . . .}), . . , 2mo + 1}). 



3.0.2 Transfer matrix 

Infinite volume. For z G C\{0}, consider the equation (U — z)ip = (dropping the 
subscript ui in the notation). In what follows we always consider z ^ 0; an analysis for 
z = can be done in a similar way, but is not of any use in this work. Due to the band 
structure of U, we obtain a recurrence relation for the components of vectors solving the 
above equation. Denote by tpn the component of tp along the basis element e„. A vector ip 
solves {U — z)tjj = if and only if its components satisfy 



V'2n+l 

where the transfer matrix is given by 



V'2n-1 
i>2n-2 



Tz{'^2ni^2n-l) — 



Zt 



We have 



detT^ 



^2gi(u>2n-l+W2n) + j.2 

-rt 



■,-i("2n-'^2n-l) 



(3.3) 



(3.4) 



(3.5) 



so equation (3.3) can be inverted. It follows that if we choose V'l and ■00 then all components 
of V' are fixed, by (3.3). A priori therefore, each eigenvalue could be doubly degenerate; 
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but this docs not happen as wc sec from the following argument. Suppose that ■i/' and x 
are two eigenvectors with the same eigenvalue z. Then we have 



V'2n+l 
V'2n 



X2ra+1 
Xin 



V'o 



Xi 
Xo 



(3.6) 



where 

T^{uj,v) := Tz{u)2n,l^2n-l) ■ ■ ■Tz{uJ2,i^l)- (3.7) 

The determinant of the matrix on the left side of (3.6) tends to zero as n — >■ oo, since tp 
and X are eigenvectors and hence belong to l^. On the other hand, \detTz{oj,n)\ = 1 for 
all n. We conclude that the columns of the matrix to the right in (3.6) are multiples of 

each other, and hence so are the vectors tp and x- 

Solutions (f^. Take z £ C from the resolvent set of U. Then there are unique (up 
to multiplication by a scalar) solutions </3^ G to {U — z)ijj = 0. Here, ti^ — {tp & l'^ : 
^^yi IV'inP < co}. To see existence we set ipo = [U — z)^^eo, so that {U — z)'ipo = cq, 
and we choose (f^ = [ipo]k for k > I, ip'/^ = [ipo]k for fc < — 1 . We then define (p^ for the 
remaining components k by applying the transfer matrix. Uniqueness of (p^ is shown just 
as in the above argument following (3.6). 
Similarly to (3.6), we have 



'fitn 
"ftn-l 



V2n- 



= T{oj, z, n) 



ft 



^1 



(3.8) 



for a matrix satisfying | det T(w, z, n)| = 1. The columns of the matrix to the right are 
linearly independent, for otherwise y+ and would be multiples of each other. This in 
turn would mean that p^ G P ^ which is not the case since z is not an eigenvalue of U. By 
taking the determinant in (3.8) wc see that for all n € Z, 



Semifinite volume. Consider U'^ acting on /^({2no, . . .}). 
z)xjj = recursively and obtain for n > hq + 2 



(3.9) 

We solve the equation {U~^ — 



tp2n 



= T+{oj,n) 



^e"^2„o+i(^e'"=^" 




r) 



V'2no+2 
V'2no 



and 



t 



-{Z 



re 



2no 



(3.10) 



(3.11) 



Here, we have set T+{u!,n) = Tz{i^2n,'^2n-i) ■ ■ ■Tzii^2no+4,i^2no+3)- A similar argument 
as in the infinite volume case shows that each eigenvalue has multiplicity one. 

Lyapunov exponent. Products of transfer matrices yield generalized eigenvectors ip solu- 
tions to {U — z)4> = and their asymptotic behaviour at infinity is related to the associated 
Lyapunov exponents. The following properties of the Lyapunov exponent will, in turn, be 
needed in the proof of Lemma 3.8 below. 

Theorem 3.1 Let /i he absolutely continuous with density r € L°°(T) and having support 
with nonempty interior. Then, there is an e > such that for every nonzero z € C with 
1 — e < |2;| < 1 + e, the limit 



lim -l0g\\T^{uj2n,i^2n-l) ' 
n— >oo n 



■T^{uj2,i^l)\\ = liz) 



(3.12) 



exists almost surely, j{z) is deterministic and strictly positive. Moreover, z 7(2;) is 
continuous. 
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Remark: The limit 7(2) is called the Lyapunov exponent. 
The proof of this theorem is given in appendix. 



3.0.3 Green's function 

Infinite volume. Let z 7^ be in the resolvent set of U, and let ip^ be the unique solutions 
of {U — z)(fi^ = in Zj.. We define the resolvent by 

= (f/„ - z)-' (3.13) 



(leaving sometimes out the superscript ui). Using the explicit form of U we solve the 

equation [{U — z)Gzen]k = ^k.n (Kronccker symbol) for Gze„ and hence obtain the matrix 
elements of the resolvent. The following is the result. For all n <E Z, we have 



and 

{eu,G-e2n-.) = t^f-' (3.15) 

Notice that these matrix elements are well defined due to (3.9). 



Semifinite volume. One can easily construct a vector with components (^^, fc e Z, 
satisfying [U — z)(p'^ = and [(U~^ — z)ip'^]k = for all A; > a := 2no. Here, is acting 
on /^({2no, •.•}). The components of are obtained by the transfer matrix for a suitable 
'initial condition' ip^^^ and y^^^+j. 

In the same way, one can construct ip^ with components ip\, k ^ T,, satisfying (U—z)(p^ = 
and [{U~ — z)ip'']k = for aU A: < 6 := 2mo. Here, U~ acts on l'^{{. . . , 2mo}). 

The Green's functions on the semi-finite spaces are defined by 

Gt = {U^-z)-\ (3.16) 

for \z\ 7^ 1 and where the operator is acting on the P space over the appropriate half-line 
of integers. Proceeding as in the infinite-volume case, one shows that (efe,G+e2n) and 
(efe,G+e2n-i), for k > 2no, are given by the r.h.s. of (3.14) and (3.15) respectively, where 
(p~ is replaced by (^". Furthermore, (e/c, G~e2n) and (cfc, Gje2n-i), for k < 2mo, are given 
by the r.h.s. of (3.14) and (3.15) respectively, where (p~^ is replaced by (p^. 

Finite volume. On Z^({2no, . . . , 2mo}) the reduced unitary is given by (3.2). The matrix 
elements of the Green's function (;;P"0'2™o] _ ^jj[2no,2m.o] _ ^-j-i gagj^y obtained in terms 
of the vectors (p"', <p^ defined in the previous paragraph. We mention explicitly only the 
formula 

(e2no,GP"'"^™''le2™o) = ~~ ^JT-^, -^^Lo-i^2„o- (3-17) 

^ ' ^2mo-1^2mo ^2mo-1^2mo 



3.1 Bound on fractional moments of Green's function 

For I -2 1 7^ 1 denote the matrix elements of the resolvent of f7u> by 

{ek,{U^-z)-^ei) = G'^,{k,l)- (3.18) 
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Theorem 3.2 Assume that the random variables ujk, k 1,, are iid, with distribution 
d/i(a;o) = T{(jJo)duio with r € L°°{T). Let < s < 1 and < e < 1. Then there exists a 
constant < C(s, e) < oo such that 

E[|G^(fc,Or] <C(s,e), (3.19) 

for all ze{CeC : 1 - e < |C| < oo, |C| 1} and all k,leZ. 

Remark: The result holds for the semifinite volume restrictions Gf{k,l) as well. 
The proof of this result is given in Hamza's thesis (Theorem 6.1) and [16]. It only depends 
on the fact that is the product D^S, where S is unitary and = diag(. . . , Wfe, w/j+i , . . .) 
has the specified random properties. 

3.2 Exponential decay of fractional moments of Green's function 

Theorem 3.3 Assume that the random variables uJk, k G Z, are iid with distribution 
dfi{u}o) = r(wo)dwo with r S L°°(T) and that the support of fi contains a non-empty open 
set. Let < e < 1. There exists an s with < s < 1/2, and there exist constants 
< C < 00, a > 0, such that 

E[|G^(fc,Or] < Ce~"l'^-'l, (3.20) 
for all z £ C satisfying \z\ ^ 1 and < |^;| < 1 + e, and for all k,l & Z. 
The proof is done in three steps that we now describe. 

3.2.1 Reduction to even matrix elements 

Theorem 3.4 Suppose that \k — l\ > 4 and let m,n be the unique integers s.t. k e 
{2m, 2m + 1} and I e {2n — 1, 2n}. Then we have 

|G^(2m-2pi,2n + 2p2)|. (3.21) 

pi,P2e{o,i} 
It follows that for any s > 

\G':{k,l)\' <C{z,ty J2 |G^(2m-2pi,2n + 2p2)r, (3.22) 

pi,P2e{o,i} 

where the constant C is that in front of the sum in (3.21) and where we have we used that 
(a + by < [a-' + ¥] for a, 6 > and < s < 1. 

Proof of Theorem 3.4- We write simply G for in this proof. We first show that for 
k ^ n, 

G(2A: + 1, 2n) = -e''^^''-'G(2k - 2, 2n) ^—G(2k, 2n), (3.23) 

r rt 

and that for A; > 2n + 2 or A; < 2n — 1, 

G(fc, 2n - 1) = -e+''^^''-iG(fc, 2n) - -^^G{k, 2n + 2), (3.24) 



|G^(fc,OI< 



\z\+t 



r rt 
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where dn = [f2n-if2n ~ 'ftnV2n-i\ ^ (^-f- (3-9))- Combining (3.23) and (3.24), and using 
the fact that \-f^\ = 1, the bound (3.21) is easily obtained. The latter fact follows from 



rf„+i = - det 



det T 



^2n-l V>2n-1 



'P2n 



C^detT, 



where | dctT] = 1. 

Let us now derive (3.24) for k <2n—l. We have G(fc, 2n — 1) = ^dn'f2n'fk^ ^^'^ ^2n 
given, using the transfer matrix as in (3.3), by ip2n — {^^t^tn-i + t^y^2n-2)- ^^i^ 

relation, we replace <^2n-2 solving (3.3) (using the inverse transfer matrix), yielding the 
expression 99+ _2 = je-''^="-V^„+i + ""'Tt""' [^^e'<"^"-^+"^") + r^]y^„. It follows that 



G(fc,2n- 1) 



— !lp-i('^2n+o.2„-i)^G(A;,2n + 2) 



+ J_e-i(c.2„+C.2„-l)[_22gi(a,2n-l+C.2n) +^2]^^j^^2n- 1), 



which yields (3.24). If fc > 2n, the same arguments with 'fi2n '^^ place of (fi2n give the result. 
Expression (3.23) is obtained in an analogous manner. ■ 

3.2.2 Reduction to finite volume 

Theorem 3.5 Assume the hypotheses of Theorem 3.2, with < s < 1/2. For any pair 
of integers m ^ n let k+ = max{m, n} and k- = min{m, n} . There exists a constant 
C^{s,t,e) < 00 s.t. we have 

E[|G,(2TO,2n)|«]' < C^(s,t,e)E[|GP-'2M(2fc_,2A;+)p'*], (3.25) 

for all z satisfying \z\^\ and 1 — < \z\ < 1 + e. 

Proof. For a given m, we decompose ^^(Z) = /^({. . . , 2m — 1}) ® /^({2m, ...}), and set 
[/ = [/-©[/++ r, where U' = L><^S'^"°°'^'""^l5l"°°'^""^' (of course, the diagonal is 
restricted to the correct half-space) and where = 0^3^"^'°°^ S^"^'°°\ The operator F 
is explicitly given by 





{k,l) 


= (2m - 1, 2m - 1), (2m - 1, 2m) 




{k,l) 


= (2m, 2m- l),(2m,2m) 


-te-'"^^"^-\ 


{k,l) 


= (2m - 1,2m- 2) 




ik,l) 


= (2m - l,2m + 1) 




{k,l) 


= (2m, 2m - 2) 




{k,l) 


= (2m, 2to + 1) 



(efc,re;) =r(A;,Z)= { 



and r(fc, I) = Q for all other values fc, I. Denote 

:= {U- ®U+ - z)-^ = G^-°°'2'»-il ® Gi2'»>°°), 



(3.26) 



(3.27) 



then the second resolvent identity gives Gz = G™ — GzTG^. Suppose that m < n — 1. 
Then G™(2m, 2n) = Gi^'"'°°'(2m, 2n) and we have 



G^(2m, 2n) = GP™'°°)(2m, 2n) - ^ G^(2m, k)V{k, l)G'^{l, 2n). (3.28) 
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According to (3.26), the matrix elements r{k, I) vanish unless I = 2m, 2m— 1, 2to— 2, 2m+l. 
However, for / = 2m — 1, 2m — 2 we have G"^(l, 2n) = (by the block diagonal form (3.27)), 
so the only terms in the sum in (3.28) are with / = 2m, 2m + 1. It follows that (still m < n) 

\G^{2m,2n)\ (3.29) 
<[l + 2[l + {\z\+r)t-^] max |G^(2m, fc)|||GP'".°°)(2TO, 2n)|. 

L fe=2m,2m— 1 J 

To arrive at this bound we use in (3.28) the relation 



e 



(2m + 1 , 2n) = [z - re-"^^"^ ] G™ (2m, 2n) 

which follows readily from the explicit expressions of the semi-finite volume Green's func- 
tion, see (3.11). 

In a next step, we estimate |Gi^™'°°^(2m, 2n)\ in (3.29) from above by the finite- volume 
Green's function |Gi^™'^"' (2m, 2n)|. We split the space as P({2m, . . .oo}) = P{{2m,2n})® 
f{{2n -hi,... oo}), so that C/I^™."") = {/[2m,2n] ^ j7[2n+i,oo) ^ f, with 



ek,rei)=T{k,l)= < 



(r - l)e-''^2"-i, (fc, I) = {2n - 1, 2n), 
ie-'---S (A;,0 = (2n-l,2n + l), 

te-i'^2n+2^ (yt, I) = (2ti + 2, 2n) ^ ^ ^ 

-{r+ l)te-^"="+= , (A;, Z) = (2n -h 2, 2n -I- 1) 

and r(A;, I) =0 for all other values k, I. By the second resolvent identity, we have (7^'"'°°) = 
^^[2m,2„] ^ ^[2n+i,oo)^(^ _ f gP™'-)). Taking the matrix elements, with m < n, we obtain 

GP™'°°)(2to, 2n) 

= GP'".2n](2TO,2n) - ^ (e2m,G[2'".2n] GP"+i'°°)efc)r(A:,0GP'".°°)(;,2n). 



Note that due to (3.30), the sum over k is really only over = 2n — 1, 2n -|- 2, and the term 

with fc = 2n + 2 is absent since the scalar product in the sum vanishes for this value of k. 
By using additionally that G^^"''^'^\2m,2n - 1) = ze'"^"-^G^^'^'^'^\2m,2n), which follows 
from the explicit formulas for the Green's function, we obtain the bound 

|(o[2m,oo)^2m,2n)| < (l + 3|0| _max \G^^"''°°\k,2n)\\ \G^^"'^^''\2m,2n)\. (3.31) 

I k — 271)271-1-1 J 

Combining (3.31) with (3.29) yields that for m < n — 1, 

\G^{2m,2n)\< {l + 2{l + {\z\+r)t-^) max \G^{2m,k)\] 

X {l-h3|^| _max \Gf"''°°\k,2n)\}\G^^"'''^''\2m,2n)\. (3.32) 

k — 271)271-1-1 J 

We take the expectation of the inequality (3.32), use Holder's inequality and Theorem 3.2 
to arrive at the bound (m < n — 1) 

E[|G42m,2n)r]' < C^{s,t,e,z) E[|G[2'"'2"l(2m, 2n)|2«] , (3.33) 

for < s < 1/2, < e < 1, |2;| ^ 1 s.t. 1 — e < |2:| < oo, and where G^(s,e,2;) is bounded 
uniformly in compact regions of z. 
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Next wc deal with matrix elements of the resolvent with m > n+1. One can proceed as 
in the above argument, or instead use the following path. Since U is unitary, \Gz{k,l)\ = 
I (efe, {U - z)-^ei) I = I (ei, (f/-i - z)-'ek) \, and (U-^-z)-^ = -l/z-{l/z)HU-l/z)-\ 
Hence we have for m > n + 1 

|G,(2m,2n)| = j^\G,/,{2n,2m)\. (3.34) 

We want to use the boimd (3.33) on the right hand side of (3.34). The condition 1 — e < 
l/\z\ < 1 + eis equivalent to 1 — < l^j < 1 + j^. It follows that under the last 
condition on \z\ we have, for m > n + 1, 

E[\G,{2m,2n)\f < C^{s,t,e,l/z) j^E[\Gf;f"'\2n,2m)\^^]. (3.35) 

Combining (3.33) and (3.35) yields the bound (3.25), provided ^ 1, 1 — e < |^;| < 1 + e 
andl- < 1^1 <1+ ■ 



3.2.3 Exponential decay in finite volume 

We assume that the support of the measure fi contains a non-empty open set in [0, 27r). This 
implies positivity of the Lyapunov exponent, see Theorem 3.1, a result we use implicitly in 
this section. 

Theorem 3.6 There are numbers a, s satisfying a>0,0<s<l such that for all z G C, 

\z\ 7^ 0, 1 and n < m, we have 



E 



|g^2„,2m](2n,2m)|*l < Ce-"( 



m—n) 



(3.36) 



The constant C depends on s, z, t and fi, but not on m, n. It is furthermore uniform in z, 
for \z\ ^ 1 restricted to compact sets ofC \ {0} (see explicit bound in proof). 

Proof of Theorem 3. 6. The proof is based on the following two Lemmas. 

Lemma 3.7 Let < s < 1. Then there is a constant < C^(s) < oo s.t. for all \z\ ^ 0, 1, 



jg^|^[2„,2H(2n,2m)|"] < |^r"(l + |^r'')(l + 2''C^(s)) E 



<P2m 
<P2m-l 



(3.37) 



Here, ^f" is the solution satisfying the left boundary condition at a = 2n, see also before 
(3.16). 

We present a proof of this Lemma at the end of the present section. The following a 
priori bound has been adapted from the self-adjoint Anderson model situation (see [11], 
Lemma 5.1), and is given in Appendix A of [14]. The proof uses strict positivity and 
continuity of the Lyapunov exponent, which we prove in Theorem 3.1 below. 

Lemma 3.8 Let A C C be a compact set not containing the origin. There are numbers 
a > and < s < 1 (depending on A), such that the product of transfer matrices (3.4) 
satisfy, for all m> n, 

E[\\T4uJ2m,0J2m-l) ■ ■ ■T,{0J2n,CV2n-l)v\\-'] <(7e-"(™-"), (3.38) 

for any normalized vector v gC^. The constant C is independent of z G A and v. 
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A combination of these two Lemmas yields a proof of Theorem 3.6 as follows. Solving 



the equation (3.3) for ¥'2m-i gives the relation ^2to-i 
therefore 



~^[t'P^m+i+r^2m]lz, and 





1 






V2m 




. 'P2m . 


~ t 


t 




. V>2m-1 . 





<P2m 
V2m-1 

We have \\T\\ < C{1 + \z\)/t, for some C independent of z, r, t and any of the phases. Using 
this expression and definition of the transfer matrices (3.3) we obtain 



f2m 
V>2m-1 



> \\Tz{(jJ2m, <^2m-l) ' ' ' Tz{(jJ2n+2, U!2n+l)v"'\ 



c(i + N) 



(3.39) 



where we chose 



'P2n+l 
V2n 



to be the normalized vector 
1 



V = 



t 



Combining (3.39) with (3.37) and (3.38) proves the bound (3.36) and hence Theorem 3.6. 
It remains to give the 

Proof of Lemma 3.7. The components of t/?^, with h = 2m, satisfy — ^'/'Im-i + 



V2m = 0; ^'^^ before (3.16). Moreover, since </? , is defined modulo a multiplica- 
tive factor only, we set (P2rn-i = 1 (note that ^2m-i — would imply that ip^ = 0, so 
the normalization (p2m-i = 1 is possible). Therefore we obtain from (3.17) the following 
expression for Green's function: 



GP"'2H(2„,2m) = - — 



1 



(3.40) 



Note that (p^m and f2m-i cannot both vanish, since otherwise we would have = 0. It 
follows from (3.3) that the components of i^" satisfy 



r2m-2 



= Tz{U2m-2,(^2m-3) 



'P2m-3 



(3.41) 



and furthermore, that <^2m = ~ — 
on uj2m and uj, with j < 2m — 2. We have 



[i'P2m-2 - rV2m-i\ ■ Consequently, (f^^ depends only 



E 



G[2"'2™l(2n,2TO) 

2tv 



\z\-'E 



d/x(W2' 



m-1) 



ZC'"^'"-^'^2m-ll' 



(3.42) 



where E is the expectation over all coj, j = 2n, . . . , 2m — 2 and j = 2m. The dependence 
on u)2m-i of the integrand in 



I:= dtx{0J2m-i)r-r 

Jo W2r, 



is concentrated exclusively in e'"2™-i_ define the vector 



Vrr,. ■■= 



V2m 
'P2m-l 



(3.43) 



(3.44) 
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If ^2. 



then wc have / = l^sfml"^ = !km|l~^- If y^2m = then wc have I = 
= \z\~''\\vm\\~'' ■ Next suppose that both </?2m-i ^2m nonzero. We 
distinguish two cases: either \<fi2m\ — \'P2m-i\ \V2m\ < 1*^2771-11- I^^ former case we 
have \\vm\\ < + \v>2m-i\ < 2|'^2ml. and 



2m-l 
\z\-'\^^m-l\ 



2n 



(3.45) 



Here, we have used that for all < s < 1 there exists < C^(s) < oo such that for all 

/3 G C 

' "d/zH|e±--/3|-^<C7^(s), 



/ 

Jo 



see [17]. Next we consider the case \<fi2m\ < Ivim-ili which case we have \\vm\\ < 

k2ml + k2m-ll<2k2m-l|,and 



\V>2m-l\ 



« / d^O 

Jo 



W2m-lj 



^ — i-,^0, I.JO, 

^ ^2m/'P2m-l\ 



Combining these estimates, we see that in any event, 

I<{l + \z\-'){l + 2'C^{s))\\Vm\\-'. 



(3.46) 



(3.47) 



This bound, together with (3.42), yields (3.37) with E instead of E. But both expressions 
are the same since docs not depend on L02rn-i- 

This completes the proof of Lemma 3.7 and with that the proof of Theorem 3.6. ■ 



3.3 Proof of Theorem 3.3 

Let s, a, C bo as in Theorem 3.6. Combining the latter theorem with Theorem 3.5, we 
obtain the bound 



|G,(2n,2m)|^/2 



< Cic 



— 2a\m — n\ 



(3.48) 



for all l^l ^0,1 satisfying the bound indicated in Theorem 3.5. The constant Ci and all 
further constants Cj introduced in this proof depend on z, s, e and t and are uniform in z 
in compacts of C \ {0}. We use the inequality (3.22) after Theorem 3.4 to arrive at the 
following bound for |A: — Z| > 4, 



E 



\G,{k,l)\'/^]<C2 E[|G^(2m-2pi,2n + 2p2)|'/' 

Pi,P2e{o,i} 



(3.49) 



Next, since \m — pi — {n +P2)| ^ l^i — n\ — 2, and \k — 1\ < \2ni — 2ti| + 2, combining (3.48) 
and (3.49) gives 



E 



\G,{kJ)\ 



sl2 



<Cse 



-a\k-l\ 



(3.50) 



provided |fc - /| > 4 and \z\ ^ 1, < |2;| < 1 + e. Finally, if |A; - /| < 4, the bound (3.50) 
is implied by Theorem 3.2. This completes the proof of Theorem 3.3. ■ 
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3.4 Proof of Theorem 2.2 



Exponential decay of fractional moments of Green's function implies dynamical localization 
for band matrices of the type — D^S, as shown in the following result. 

Theorem 3.9 ([16], Theorem 3.2) Assume that the random variables uJk satisfy the 
conditions of Theorem 2.2 and that for some s € (0, 1), C < oo, a > and e > 0, 

E(|G^(A:,/)|'') < Ce""!''-'! 

for all k,l G Z and all z gC s.t. 1 — e < |^| < 1. Then there is a C < oo such that 



E 



sup \{ek,f{Uu,)ei) 
/eC'(S),||/|U<i 



for all k,l G Z. 



With the identification of the basis elements e2k = | t) <8i e2k+i = \ i) ^ see after 
(2.21), we immediately obtain that Theorems 3.3 and 3.9 imply Theorem 2.2. ■ 



4 Temporal disorder 

For the sake of comparison, we briefly consider in this section the case where the disorder is 

introduced in the model through the time variable. This means that at each time step, the 
evolution of the coin variable is randomly chosen within a set of iid unitary 2x2 matrices 
{C'fejfeez- Therefore, the random evolution after n steps reads 

U{n,0) = UnUn-i---U2Ui, where Uk^S{Ck<»l). (4.1) 

As we will see, our choice of random unitary matrices {Ck}kez 



Ck. 



e-'"it -e-'"ir 



(4.2) 



e~'"fc r e~"^h: t 
with {Wfc }feez U {w^jfeez iid subjected to the condition 

E(e''^^) =E(e'"^") =0, (4.3) 

naturally leads to the study of (classical) persistent random walks. For notational reasons 
which will be clear below, we change notations tot = +l.i = — Iso that 

It) = 1 + 1), i;> = i-i>, Pt = p+i, Pi = p-i- (4.4) 

We first state a deterministic result dealing with expectation values of position opera- 
tors at time n. 

Let X denote the position operator on (g) i^(Z) defined by {Xil)){x) = (I (g) x)ip{x), 
X Gil, on its maximal domain 

r>= .j V'S C2 0/2(Z), s.t. ^ ^ \\{P^®x)i^{x)\\l^ <Oo) . 

xezcre{+i,-i} 
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For / : Z C, wc define the operator F{X) on ® by iFip){x) = (I (g) f{x))'ip{x) 

on its maximal domain Dp via the spectral Theorem. For any S (g such that 
{/(n, 0)^0 € -Df, we note 

(vfo, U{n, OrF{X)U{n, 0) Vo) = (i^(^))^„ (4.5) 
Explicit computations with P±i — \ ± 1)(±1| yield the following Lemma: 
Lemma 4.1 With the notations above, 

n 

n 

= Y Jk{n)^\x + k){x\ (4.6) 

xG^ k=—n 

where 

Jk{n)= Y P.nCn---Pa,CieM2{C) (4.7) 

satisfies Jk{n) = if k and n have different parities or if k > n or k < —n. 
Moreover, if = ® \0) , we have for any / : Z i-)- C and all n € N, 

n 

{F{X))^M= ^ /(/c)(v5o,Jfe*(n)Jfe(n)v5o)c2, (4.8) 

k=—n 

where Wk{n) = {^o,J^{n)Jk{n)(pQ)^-2 satisfies 

n 

Wk{n) >0 and ^ Wk{n) = \\<fio\\^. (4.9) 

k——n 

Remarks: 

i) If the initial coin vector ipo is normalized, which we assume from now on, then the 

quantum mechanical expectation value of the operator F(X) at time n coincides with 
the expectation value of the function / with respect to the classical discrete probability 
distribution {H^fc(?i)}fce{-Ti,...,Ti} on {— n, ...,n} C Z. The quantity Wfc(n) is interpreted 

as the probability to reach site fc € Z in n steps. 

ii) The probabilities {Wk{n)}ke{-n,...,n} are actually </Jo-dependent random variables for 
randomly chosen coin operators {Ck}kei. with arbitrary distribution. Taking expectation 
with respect to the distribution of the Cfc's, we get a new discrete probability distribution 
on {-n, . . . ,n} C Z given by {wk{n)}ke{-n,...,n} with 

Wk{n) =E{Wk{n)), k G {-n, . . . ,n}, V n e N, (4.10) 

with same interpretation in terms of a classical random walk on Z. 

We shall focus on the expectation value of the quantum mechanical moments of the 
position operator for our choice (4.2) of {Ck}k& under condition (4.3), for a normalized 
initial condition of the form tpo = ipo^\0), i.e. on 

n 

E{{X^)^„{n)) = Y k^Wk{n) as n^oo. (4.11) 

k=—n 
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We remind the reader that a persistent or correlated random walk on Z is determined 
by a probabihty p that the walker moves one unit in the same direction as in the previous 
step, and a probability 1 — p that the walker moves one unit in the opposite direction to 
the last step. See e.g. [28]. 

We let Sn = J2^=i '^3 with aj G {—1, +1} be the random walk defined by P(-S'„ = k) = 
Wk{n). This random walk is characterized as follows. 

Proposition 4.2 Assume the matrices {Cfejfegz are given by (4.-2) and (4-3) holds. Take 
V'o = Vo <8i |0) with ipQ = a\+) + 13\—) G normalized. Then Sn is a persistent random 
walk with parameters 

P((Ti = +1) = lapt^ + - 23fJ(a/3)rt := a 

P(cri = -1) = japr^ + \p\^t^ + 2^{aP)rt := b = 1 - a 

F{aj = +l|crj_i ^ +1) = P(fTj = -l\aj-i = -1) = t"^ 

¥{aj = -l\aj-i = +1) = Piaj = +l\aj-i = -1) = = 1 - t'^, 

for j > 2. 

Proof: For n € N, let ct = (cti, (72, . . . , (t„) € {—1, +1}". We need to consider 

Wk{n)= Yl E{{ipo,C*,P,,C^P,>---C*P,,^P,„Cn---Pa,Ci^o)), (4.12) 

o-.tr'ef-i.+i}" 
E"=i<'i=Ei=i<'5='= 

where cr„ = a'^. With P„ = \a){a\, a = ±1, and reorganizing the product, the scalar 
product under the sum equals 

n 

where the first two factors can be further expanded as 

(a + ^ K,Ci-)) (a(ai , Ci +) + , Ci -)) . 

By independence, the expectation factorizes and, with the notation {(j,Ct) = Ccr,T-, cr,T G 
{+1, —1}, we immediately get from (4.2) and (4.3) 

lE(a^a,,a,_J = if a'j j^aj, V j > 2. 

This, together with the conditions cr„ = a'^ and J^j^i — Sj=i "'j = k in (4.12), imposes 
aj = a'j for all j > 1. Hence, (4.12) reduces to 

Wk{n) = (4.13) 

n 

^ E(|a|2|a„+i|2 + + 25R(a^a„_ia::;T)) n E(|a,..,_i I'), 

<^e{-i,+i}" j=2 

where 

E(|C+i,+ip)=E(|C_i,_ip)=t2 
E(|C+i,- i|^)=E (|C_i,+i|^) = 

E(C_i,_iC_i,+i) = -E(C+i,_iC+i,+i) = rt. (4.14) 
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The result then follows with the definition of Sn = X^"=i (^j ■ ^ 

Persistent random walks or correlated random walks are well known and have been 
studied in many details and greater generality, ours being the simplest instance. See e.g. 
[28], [12] and the references therein. In particular, when r = t = l/-\/2, the persistent and 
symmetric random walks are equivalent. Also it is known that the first moment is finite 
and that the second moment is proportional to n for n large. This leads to the following 
diffusive behaviour: 



lim 

n^oo 



(4.15) 



For completeness, we provide below a simple proof of the fact that all moments of the 
persistent random walk display a diffusive behaviour, a statement that we couldn't find as 
such in the literature, although certainly well known. 

Proof of Proposition 2.4. We assume the hypotheses of Proposition 4.2 and use the 
familiar setup of generating functions together with a classical scaling argument. 
We consider only the situation r which differs from the usual symmetric random walk. 

Let ti'^(n) be the conditional probabilities 



w^{n) =¥{Sr,^k\an = ±l) s.t. 
Thus we have, for n > 1 and Ifcl < n 



wl{n)+w^, (n) = Wk{n). 



w^{n + 1) 



Wk{n+l) = i^w^^^ (n) + r^w^+i (n) 



(4.16) 



with 



wt{l) = a, u;f(l) = 6. (4.17) 
Moreover, w^(n) = if |fc| > n. We introduce the generating functions and ^„ by 



k=—n 

*„(^) = ^+{z) + ^-{z), yzec. 

As a consequence of (4.16), introducing ^n{z) = ($^(-2), $"(2;))^, we have 

^n+l{z) = M{z)^„{z), with 



(4.18) 
(4.19) 



M{z) 



and $i(-2) 



he- 



This allows to determine explicitly $„(z) and 



M"-i(2)$i(2)), n>l 



(4.20) 



and, in turn, all moments of the probability distribution {wfc(n)}fegz- 



We consider now the diffusive scaling introducing the macroscopic time variable N = 
nr, where r >> 1, r e N and the macroscopic space variable K = ^/rk, such that 
K/^/N — kj^fn remains finite. Expecting a probability distribution {uik{p^"\keTL asymp- 
totically invariant under this scaling, we are led to the study of the generating function at 
z = yl\fT in the limit r ^ 00: 
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Lemma 4.3 



k=-Tn 



(4.21) 



uniformly in y in any compact set of C 



Since the functions of y involved are entire and the convergence is uniform in compact sets, 
wo can differentiate the above identity w.r.t. y as many times as we wish. In particular, 
differentiating L times with n = 1 and setting y = immediately yields the following 
Corollary which ends the proof of the Proposition 2.4: 



Corollary 4.4 



lim 



^ [ikY 



L/2 



(4.22) 



where Hl denotes the Hermite polynomial Hl{z) = {—l)^e^^/^ (^) e 
Proof of Lemma 4.3. Note that the left side of (4.21) equals 



lim *Tn(j//v^) = lim 



.M-"-i(y/v^)$i(t//v^) 



(4.23) 



so that we are interested in small values of 2: = y/y/r. For l^;] small enough, we compute 
the spectrum of the analytic matrix M{z) 



(t{M{z)) = {X,{z),X2{z)} 



with distinct eigenvalues 



Ai (z) = t"^ cos(z) ± ^/t^ cos2(z) + l-2t^. 
Hence, for \z\ small, there exists an invertible matrix R{z), analytic in z, such that 



M"(2) = R-\z) 
For \z\ small enough, we have 



A^(z) 
X^iiz) 



R{z). 



(4.24) 

(4.25) 

(4.26) 



t' 



\,{z) = l-z'— + Oiz') 



X^{z) = {e-r'')+z 

whereas 

M-i(O) = M-i(O)* 

and 



,tHe-r^ 

2r2 



+ 0{z^), 



(4.27) 



1 



2 _^2 

r2 ^2 



satisfies M-^(O) 



" 1 ■ 




■ 1 " 


1 




1 



R{Q) = R*{Q)=R-\Q) = — 
Since (4.27) yields for any y € C and r large enough 



1 1 
1 -1 



Ai(y/\/?)"" = e-"'''^+0(nyVr) 



if _ r-2)^"e"f'^ + 0{ny^/T), 



(4.28) 
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we eventually obtain, since — < 1 and a + b= 1, 



lim 



rn—l 



(2//v^)$i(|//v^) 



-ny 





■ 1 ■ 




' a + b ' 









5 


a — b 


> 



(4.29) 



A Appendix 

A.l Proof of Lemma 2.1 

Introducing the discrete Fourier transform T : i^([0, 27r), C^) (g) by 



(A.l) 



we get that U is unitarily equivalent to a multiplication operator by a matrix V{x). 



J^-'^UF = multy(a;) = mult 







C. 



(A.2) 



The matrix V{x) is entire in x, and unitary for x real, which means that its eigenvalues 
{a\{x),a2{x)} and eigenprojectors {Pi(a;), ^2(2;)} are also analytic functions of a; S M, 
even at the possible crossings of eigenvalues for a; e M. Moreover, for any n e Z, 



V^{x) = Pi{x)a'l{x) + P2{x)a^{x). 



(A.3) 



Similarly, the position operator K — (I<S> k) is unitarily equivalent to differentiation w.r.t. 
X (on its natural domain): 

T-^KT = -id^. (A.4) 

In particular, for e f{Z) such that J""!* = / S L^qq, 27r), C^) and * in the 
domain of = (I (g) /c)^, we have for all n e Z, 

/■27r 

{K^Mn):={^,U-^K^U^^)= \\d4V^{x)f{x))\\l.dx. (A.5) 

Jo 

Now, it is not difficult to see by explicit computations that this quantity behaves as for 
n ^ 00, unless the analytic eigenvalues {ai{x),a2{x)} of V{x) are independent of x. We 
have 

ai{x)a2{x) = det V{x) = det C and ai{x) + a2{x) = tr V{x) = e~"a + e"d. (A.6) 
Hence the eigenvalues aj{x) are independent ofa;iffa = rf = 0. H 

A.2 Posit ivity and continuity of Lyapunov exponent 

This Section is devoted to the proof of Theorem 3.1. 

Remember that wc take /i to be absolutely continuous with density t <E L°°{T) having 
support with nonempty interior and that the transfer matrices T;^{9,r]) are given in (3.4) 
(with u)2n = and u)2n-i = v)- 
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To prove Theorem 3.1, we use Furstenbcrg's theorem, which is developed for real square 
matrices. We thus map Tz € M2(C) (square 2x2 matrices with complex entries) into 
t{Tz) e M4(K) using the bijection 



where 



a b 




c d 








" 1 


I = 






IRea + JIma 
/Rec + JImc 



/Re6 + Jlmb 
IRed + JImd 



(A.7) 



and J = 







We refer to [8] for more detail on this transformation. In particular, we have 

\\TiA)\\ = V2\\Al 



(A.8) 



where the norms arc given by HXjp = Tr(X*X), for X € M2(C) or X G M4(M), and 
that if A e M2(C) satisfies |det(A)| = 1, then |detT(A)| = 1. Due to (3.5), we have 
I det TiTzie,rj))\ = 1, for all ^ ^ and all 0,7]. 

Relation (A.8) together with the fact that t{AB) ~ t{A)t{B) shows that the state- 
ment of Theorem 3.1 is equivalent to lim„^.oo ^\\T{Tz{9n,rin)) ■ ■ ■ T{Tz{9i,r]i))\\ = 7 almost 
surely, for the same deterministic 7 > as in Theorem 3.1. 



(A.9) 
) with 



The measure on T induces a measure on M4(1R), supported on the subset 

M ■.= {T{T,{0,r])),e,r]esnpp]i}. 

We call the induced measure again fj.. Let Fi, 12, • • • be iid random matrices in IV 
common distribution fx. If the integrability condition 

E[max{log||Fi||,0}] < 00 

is satisfied, then the upper Lyapunov exponent 7 € R U {—00} is defined as 



(A.IO) 



7:= lim -E[log||F„---yi||]. 

The theorem of Furstenberg and Kesten ([7], Theorem 4.1) states that if in addition the 
matrices Yj are invertible {Yj e GL4(]R)), then 



lim -log||F„ 

n— )-oo 77, 



.•yi||=7 



(A.11) 



almost surely. In our case, Yj = T{Tz{Oj,rjj)) is invertible, and ||Yi|| < C is uniformly 
boimdcd in 0i,?7i (see (3.4)), so that (A.IO) is trivially satisfied. This implies that (3.12) 
holds almost surely, with a deterministic 7. Moreover, since S = Tz{0n, Vn) • • • Tz{0i,r]i) is a 
2x2 invertible matrix, we have US'!! = HS"""^!! > 1, and therefore 7 > 0. The remaining part 
of the proof of Theorem 3.1 consists in proving that 7 is strictly positive and continuous. 

Let Qfj^ C GL4(IR) be the (multiplicative) group of matrices generated by the transfer 
matrices Tz{9, r)), where 0, rj vary throughout the support of the measure ii. Here, z is fixed 
and not displayed in Q^. 

Theorem A. 1 (Furstenberg, [7] Thm. 6.3) Suppose that is strongly irreducible and 

non-compact. Then the upper Lyapunov exponent associated with any sequence of random 
matrices Yi,Y2,. . . in SL4(M) n supp/x, iid with common distribution fi, is strictly positive. 
This means that (A. 11) holds with 7 > 0. 
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Now we show that is strongly irreducible and non-compact. 
Lemma A. 2 /f supp/i contains two distinct points, then is not compact for all z ^ 0. 
Proof. For fixed z = -Re'", we have 



~R 



-rt 



and a direct calculation shows that 

T,{n,e)T,{T},n)-' = {T,{e,e)-'T,{e,r,)r 







:(ei(«-'7) 



Both operators Tz{ri,6)Tz{ri,r]) ^ and Tz{6,6) ^Tz{6,r]) belong to Q^, and hence so does 
their positive-definite product, 



< M := T,{r],0)T,{T],ri)''T,{e,e)-'T,ie,r]) = 1 - 







Note that the r.h.s. docs not depend on z at all. Since the determinant of the last matrix 
is — frje'^^"''-' — <0,Gij, contains a non-negative matrix M having an eigenvalue strictly 
larger than one. It follows that ||M"|| — oo as n — )• oo. Thus the group generated by the 
T{Tz{9,r])), 0,r] G supp/x, is not compact. ■ 

Finally we prove strong irreducibility of CJ^. Let J be a non-empty open interval, 
J C supp/U. Since the subset of matrices M' obtained from M, (A. 9), by restricting 
^, 77 e J, is a continuous image of the connected set J x J c M^, A4' is a connected set 
in M4(R). Strong irreducibility of Ai' (which implies strong irreducibility of A^) is then 
equivalent to irreducibility of A4' , see [7] Exercise IV. 2. 9. 

Lemma A. 3 The only subspaces V CM.'^ invariant under the action of M' are V = {0} 
and V = M.'^. Hence A4' is strongly irreducible. Since A4' C G^, we have that is strongly 
irreducible. 

Proof. We just need to show that V = {0} and 1/ = are the only invariant subspaces 
of M.'. Let z = i?e'" ^ be fixed. One easily finds that 



TiTz) = - [Ti cos ip + T2smip + T3 cos x + Ti sin x] 



(A.12) 



where (p = a + rj, x = ct + O vary in the open interval a + J, and the matrices Tj are given 

by 



10 

10 







To = R-' 










-rt 















-rt 


-rt 





t^ 










-rt 





t^ 






1 










-1 

























1 



















-r2 







rt 









-rt 








rt 







-f 


-rt 







t^ 






By taking derivatives in the angles, wc see that if V is invariant under M' , then V is 
invariant also under the action of — Ti sin <p -|- T2 cos </? (and — T3 sin x + cos x) , and hence 
(by again differentiating) V is as well invariant under the action of T\ cos </? -|- T2 sin (and 
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Ts cos X + ^4 sin x) . Therefore, — Ti sin^ ^ + T2 sin ip cos and Ti cos^ + 72 sin 9? cos (p leave 
y invariant, and hence so does T\ (the difference) and consequently T2 as well. Similarly 
one sees that T^^Ti leave V invariant, too. Hence any subspace V invariant under M.' 
must be invariant separately under Ti, T2, T3 and T4. 

Consider first V with diml^ = 1, i.e., V = {v) (real span). It is easy to see that the 
only possibility for v resulting in a F invariant under Ti + T2 and Ti — T2 is = ae^ + ^64 
(canonical basis elements of R^). But cither such V is not left invariant under the action 
of T3. Consequently, no one-dimensional subspace is invariant under A^. 

Consider next V with dimV = 2. Since T3 is real symmetric and must leave V invari- 
ant, V must be spanned by two eigenvectors of T3. One easily finds that T3 has eigenvalues 
{0, 1}, both twice degenerate, that a basis for the kernel of T3 is {[1, 0, r/t, 0]*, [0, 1, 0, r/i]*} 
(transpose), and a basis of the eigenspace with eigenvalue 1 is {[1, 0, —t/r, 0]*, [0, 1, 0, —t/r]*}. 
There are three possible cases: 1. both eigenvectors spanning V belong to the kernel of T3, 
2. both eigenvectors spanning V belong to the eigenspace with eigenvalue 1 of T3, or 3. 
one eigenvector belongs to the kernel of T3 and the other one belongs to the other spectral 
subspace. Either of these three cases can be analyzed separately, and one finds that none 
of the thus formed spaces V with dim V = 2 is invariant under all of the Tj, j = 1, 2, 3, 4. 
In conclusion, no two-dimensional subspace is invariant under A4. 

Consider now V with dim = 3. Then V^-*" has dimension 1 and is invariant under T3, 
since the latter is real symmetric. Hence V-^ is spanned by one of the eigenvectors of T3. 
In the same way, one sees that V-^ must also be invariant under Ti, and one finds easily 
that this implies that V-^ = {0}, a contradiction to dimF = 3. This shows that there is 
no three-dimensional subspace invariant under A4. Lemma A. 3 follows. ■ 

This proves all assertions of Theorem 3.1, except for the continuity of 2; 7(2). How- 
ever, the latter has been shown to hold in Section VII of [15] (Theorem 7.1). ■ 
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